home *** CD-ROM | disk | FTP | other *** search
/ Space & Astronomy / Space and Astronomy (October 1993).iso / mac / TEXT / SPACEDIG / V10_0 / V10_089.TXT < prev    next >
Internet Message Format  |  1991-07-08  |  19KB
  1. Return-path: <ota+space.mail-errors@andrew.cmu.edu>
  2. X-Andrew-Authenticated-as: 7997;andrew.cmu.edu;Ted Anderson
  3. Received: from beak.andrew.cmu.edu via trymail for +dist+/afs/andrew.cmu.edu/usr1/ota/space/space.dl@andrew.cmu.edu (->+dist+/afs/andrew.cmu.edu/usr1/ota/space/space.dl) (->ota+space.digests)
  4.           ID </afs/andrew.cmu.edu/usr1/ota/Mailbox/4Z8nFau00VcJ4CUU5F>;
  5.           Fri, 29 Sep 89 05:22:47 -0400 (EDT)
  6. Message-ID: <kZ8nFEm00VcJICSk4x@andrew.cmu.edu>
  7. Reply-To: space+@Andrew.CMU.EDU
  8. From: space-request+@Andrew.CMU.EDU
  9. To: space+@Andrew.CMU.EDU
  10. Date: Fri, 29 Sep 89 05:22:25 -0400 (EDT)
  11. Subject: SPACE Digest V10 #89
  12.  
  13. SPACE Digest                                      Volume 10 : Issue 89
  14.  
  15. Today's Topics:
  16.          Frequently Asked Questions, Higgins version
  17.        Re: NASA Headline News for 09/20/89 (Forwarded)
  18.          Re: ...hospital radiation machine...
  19.    Erich von Daniken (was Re: Analysis of Martian "Face" Announced)
  20.                 NEW SUBSCRIBER
  21.         C orbital prediction subroutines, and a query
  22. ----------------------------------------------------------------------
  23.  
  24. Date:     Thu, 21 Sep 89 13:11 CDT
  25. From: Bill Higgins-- Beam Jockey <HIGGINS%FNAL.BITNET@UICVM.uic.edu>
  26. Subject:  Frequently Asked Questions, Higgins version
  27. Original_To:  Orig_To! SPACE,SPENCER,EUGENE@AMELIA.NAS.NASA.GOV
  28.  
  29. I have been irritated every time Eugene posts something about Frequently Asked
  30. Questions, for two reasons:  A lot of F.A.Q.'s don't seem to be on his list,
  31. and some of the things he's answering don't seem to be Asked all that
  32. Frequently.  At last, after reviewing a few months' worth of traffic (and
  33. reading the Digest for a few years), I drew up my own list of F.A.Q.'s.  You
  34. are free to consider this just as eccentric as Eugene's.   Maybe we can kick
  35. this around a bit...
  36.  
  37. I notice the absence of questions about the Soviet programs.  I guess these get
  38. dealt with pretty well by such people as Henry Spencer and Glenn Chapman as
  39. they come up; I don't feel that there are a lot of recurring ones.
  40.  
  41. "HOW DO I GET'S"
  42.  
  43. [Canned answers to these would be good to have around, and maybe reposting them
  44. infrequently would be helpful. Note that many good definitive answers to them
  45. have been posted in the recent or distant past. Maybe we should collect such
  46. answers someplace.  Mike Smithwick once did a good job on NASA Select, Eugene
  47. wrote about astronaut selection, Ken Jenks catalogued aerospace employers, etc.
  48. Some answers are just a bit of library work.]
  49.  
  50. --How do I get access to the International Astronomical Union bulletins?
  51. --Where can I get free material about NASA/ESA/other space programs?
  52. --How do I apply to be an astronaut? How do I get 'em to select me?
  53. --Will there be any launches while I'm visiting Florida?
  54. --Where is the best place to watch a Shuttle launch from?
  55. --Tell me about NASA Select, the satellite channel.
  56. --How do I get Voyager/Viking/etc. images in machine-readable form?
  57. --How do I receive and interpret weather satellite pictures?
  58. --Where can I get software to simulate the entire Solar System?
  59. --I am writing software to simulate the entire Solar System.  Where can I get
  60. positions and velocities of all the planets and satellites?
  61. --How do I contact L5/NSS/SSI/PlanSoc/World Space Foundation/SEDS et al?
  62. --Is there an astronomy newsgroup on Internet/Bitnet/etc.?
  63.  
  64. PRETTY SPECIFIC QUESTIONS
  65.  
  66. [These keep recurring.]
  67.  
  68. --Can we just resurrect the Saturn V?
  69. --What is the orbital debris hazard?  How can it be cleaned up?
  70. --Where do Shuttle names come from?
  71. --Can Hubble Space Telescope look at the ground?
  72. --What is the current status of private launch companies? (Amroc, Space
  73. Services, Third Millenium, Pacific American, OTRAG)
  74. --Why don't those idiots at NASA put external tanks in orbit?
  75. --How would you terraform Mars? Venus?
  76. --What does a Shuttle cost?  Will private operators ever run one?
  77. --Why didn't they send Voyager to Pluto?  Or back to Earth?
  78.  
  79. "TELL ME'S"
  80.  
  81. [Some subjects keep coming up, and there is broad interest in them, rather than
  82. more narrowly focused specifics.  Perhaps our best defense is to have a list of
  83. readily available references for each subject.  Not possible in all these
  84. cases.]
  85.  
  86. --Tell me about spy satellites. (Optical, signal intelligence, radar)
  87. --Tell me about the Global Positioning System.
  88. --Tell me about amateur radio satellites.
  89. --Tell me about scramjets.
  90. --Tell me about laser launching.
  91. --Tell me about electric rockets.
  92. --Tell me about solar sails.
  93. --Tell me about electromagnetic catapults (railguns, mass drivers).
  94. --Tell me about mass extinctions, asteroid/comet bombardment, and the Nemesis
  95. (hidden companion of Sun) theory.
  96.  
  97. =============================================================================
  98. So. What have I left out?  What should we do next?
  99.  
  100.                       ______meson   Bill Higgins
  101.                    _-~
  102.      ____________-~______neutrino   Fermi National Accelerator Laboratory
  103.    -   -         ~-_
  104.  /       \          ~----- proton   Bitnet: HIGGINS@FNALB.BITNET
  105.  |       |
  106.  \       /                          SPAN/Hepnet/Physnet: 43011::HIGGINS
  107.    -   -
  108.      ~                              Internet: HIGGINS%FNAL.BITNET@UICVM.uic.edu
  109.  
  110. ------------------------------
  111.  
  112. Date: 21 Sep 89 20:13:59 GMT
  113. From: sjsca4!news@uunet.uu.net  (Greg Wageman)
  114. Subject: Re: NASA Headline News for 09/20/89 (Forwarded)
  115.  
  116. Opinions expressed are the responsibility of the author.
  117.  
  118. In article <32247@ames.arc.nasa.gov> yee@trident.arc.nasa.gov (Peter E. Yee) writes:
  119. >
  120. >This is NASA Headline News for Wednesday, September 20.....
  121. >
  122. >
  123. >As Hurricane Hugo continues to move in a northwesterly direction 
  124. >towards the U.S. East Coast, Kennedy Space Center launch 
  125. >officials have decided not to make a decision yet on moving the 
  126.                 ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  127. >Space Shuttle Atlantis back to the Vehicle Assembly Building.  
  128.  
  129.  
  130. >Meanwhile, launch officials have made a tentative decision not to 
  131.                                          ^^^^^^^^^^^^^^^^^^
  132. >demate the FltSatCom satellite from the Atlas-Centaur booster at 
  133. >the Cape.  The launch has been slipped to no earlier than Sunday 
  134. >at 4:13 A.m., Eastern time. 
  135.  
  136. It's a good thing the managers at NASA are so *decisive*.  Wow.
  137.  :-)
  138.  
  139.  
  140. Copyright 1989 Greg Wageman    DOMAIN: greg@sj.ate.slb.com
  141. Schlumberger Technologies    UUCP:   {uunet,decwrl,amdahl}!sjsca4!greg
  142. San Jose, CA 95110-1397        BIX: gwage  CIS: 74016,352  GEnie: G.WAGEMAN
  143.         Permission granted for not-for-profit reproduction only.
  144.  
  145. ------------------------------
  146.  
  147. Date: 21 Sep 89 22:48:46 GMT
  148. From: vsi1!daver!lynx!neal@apple.com  (Neal Woodall)
  149. Subject: Re: ...hospital radiation machine...
  150.  
  151. In article <Added.8Z6EOTe00UkT4D_k8g@andrew.cmu.edu> Gregory E. Gilbert writes:
  152.  
  153. >Neal Woodall mentioned in V10 #58 some unfortunate individuals who came
  154. >across an old Cobalt Therapy machine and started playing with the source.
  155. >Does anyone have a full account of this incident?
  156.  
  157. I wish i could give you more info, however it was about a year or two ago
  158. and I read a lot of newspapers and magazine.
  159.  
  160. I did happen in Brazil, several exposed people died, a few city blocks
  161. had to be decontaminated. I remember a followup story that said the bodies
  162. of the victims had to be buried in lead boxes, and the graves were checked
  163. regularly for leakage. Must have been some pretty nasty stuff.
  164.  
  165. >Bye the way, I am new to this digest what's an RTG?
  166.  
  167. Radioisotope Theromoelectric Generator. Not a nuclear reactor....the isotope
  168. of Pu used is radioactive, but not fissionable. It produces energy from its
  169. decay (heat) that is converted to electricity, but it does not support
  170. sustained chain reactions. Less power, but safer.
  171.  
  172.  
  173.  
  174.  
  175.  
  176. Neal
  177.  
  178. ------------------------------
  179.  
  180. Date: 20 Sep 89 20:17:53 GMT
  181. From: mcsun!ukc!reading!cf-cm!cybaswan!iiitsh@uunet.uu.net  (Steve Hosgood)
  182. Subject: Erich von Daniken (was Re: Analysis of Martian "Face" Announced)
  183.  
  184. In article <9189@june.cs.washington.edu> dfkling@cs.washington.edu (Dean F. Kling) writes:
  185. >    A few months ago PBS in Seattle had a British TV show that discussed 
  186. >these  kinds of ad-hoc correlations.  As an examle they took Stonehenge, 
  187. >the Salisbury Cathedral, and several local farms and drew dozens of striking
  188. >conclusions based on their geographical realtionships, sight lines, and 
  189. >long distance relationships to other English monuments, and to solstice lines, 
  190. >etc.  The obvious conclusion they drew was that the whole complex *must* 
  191. >have been part of an overall coherent design, by a single culture.   
  192. >There were just *too* many "coincidences" to be explained by mere chance!  
  193. >The lesson was that if you look hard enough for post-facto relationships, 
  194. >you will find them.
  195. >
  196.  
  197. About 15 years ago, the BBC's 'Horizon' science programme did a hatchet-job
  198. on the then-popular theories of Erich Von Daniken (author of 'Chariots of
  199. the Gods?' and other titles).
  200.  
  201. The bit of this program that sticks in my mind was where they rubbished the
  202. 'Nazca UFO Landing Strip' pictures. They did a good job of this... setting the
  203. camera up to replicate the still photo in the book. Sure, it looked just
  204. like a runway, complete with little bays for hard-standing at the sides.
  205.  
  206. Then the narrator walked into the shot. Either very small aliens built the
  207. landing strip, or the BBC have a 250 ft high narrator! The whole design was
  208. evidently no more than 10ft from end to end.
  209.  
  210. Erich Von Daniken admitted in the following interview that he didn't take
  211. the picture, but someone sent it to him and that he published it without any
  212. checks.
  213.  
  214. -----------------------------------------------+------------------------------
  215. Steve Hosgood BSc,                             | Phone (+44) 792 295213
  216. Image Processing and Systems Engineer,         | Fax   (+44) 792 295532
  217. Institute for Industrial Information Techology,| Telex 48149
  218. Innovation Centre, University of Wales,        | JANET: iiitsh@uk.ac.swan.pyr
  219. Swansea SA2 8PP                                | UUCP: ..!ukc!cybaswan!iiitsh
  220. -----------------------------------------------+------------------------------
  221.             My views are not necessarily those of my employers!
  222.            Please note the change in my username as of 1 Aug '89
  223.  
  224. ------------------------------
  225.  
  226. Date:     Thu, 21 Sep 89 10:09:30 EDT
  227. From: Michael Ryan (FSF) <mryan@PICA.ARMY.MIL>
  228. Cc: MRYAN@PICA.ARMY.MIL
  229. Subject:  NEW SUBSCRIBER
  230.  
  231. Request to be put on distribution
  232.  
  233. ------------------------------
  234.  
  235. Date: 21 Sep 89 12:33:26 GMT
  236. From: ecsvax.uncecs.edu!dukeac!tcamp@mcnc.org  (Ted A. Campbell)
  237. Subject: C orbital prediction subroutines, and a query
  238.  
  239.  
  240. Some readers of sci.astro and sci.space may have seen the Space 
  241. Flight Simulator, version 0.02, which was posted to comp.binaries.ibm.pc.  
  242. If I ever make up my mind just how to do it, I may release 
  243. the source code for the whole thing.  In the meantime, I'm going to 
  244. post the following subroutines to sci.astro and sci.space for 
  245. two reasons (a) someone might actually find them useful, and 
  246. (b) someone might want to advise me on a technical point that 
  247. is apparently beyond the abilities of this bear-of-very-little-brain.  
  248.  
  249. The following code contains three critical subroutines (in C).  
  250. or_init() sets up orbital parameters based on an apoapsis 
  251. and periapsis (in kilometers).  It presumes that the following 
  252. globals have already been set:  
  253.  
  254.     (double) fo_mass    the mass of the orbital focus
  255.                 in grams (earth = 5.98e27)
  256.  
  257.     (double) fo_radius    the radius of the orbital focus 
  258.                 in kilometers
  259.  
  260.     (double) fo_siday    the sidereal period of the orbital 
  261.                 focus in seconds 
  262.  
  263. The include file "sfsm.h" defines a structure for orbits, the elements
  264. of which should be obvious.  sfs_orbits[ n ] is an array of pointers to orbit 
  265. structures, so that sfs_orbits[ 0 ]->period would give the orbital 
  266. period in seconds for orbit 0, etc.  
  267.  
  268. The subroutine or_ssp() calculates the subsatellite point for an
  269. orbit, given the time (in seconds).  It also calculates the 
  270. current altitude.  It, in turn, calls or_kep(), which solves 
  271. Kepler's equation.  These routines may be useful for anyone 
  272. wanting to build C applications with orbital calculations.  
  273. The subroutines are based on pseudocode given in Martin Davidoff's
  274. Satellite Experimenters Handbook (published by ARRL for Radio Amateurs 
  275. using the OSCAR and other Amateur satellites), but the C implementation is
  276. my own.  
  277.  
  278. PROBLEM:  These routines seem to work well enough for orbits with 
  279. inclinations near 0 and near 180 degrees (by the way -- all angles 
  280. in these routines are in radians).  But for orbits with inclinations 
  281. near 90 degrees, there is a grotesque skewing that occurs (at the 
  282. highest/lowest points in the orbit).  As of yet, I haven't been able to figure 
  283. what causes this skewing.  
  284.  
  285. Of course, it's possible that the algorithms are correct, and that 
  286. the skewing occurs in other routines that manipulate their output in 
  287. various ways.  On the other hand, the skewing occurs so frequently 
  288. that I suspect the culprit is here somewhere.
  289.  
  290. email:  tcamp@dukeac.ac.duke.edu
  291.  
  292. /****************************************************************
  293.  
  294.     SFS_OR.C        Orbital calculation routines for SFS
  295.  
  296.             Copyright (c) 1989, Ted A. Campbell
  297.  
  298. ****************************************************************/
  299.  
  300. #include "math.h"
  301. #include "time.h"
  302. #include "vdi.h"
  303. #include "sfsm.h"
  304.  
  305. extern    double or_lfix();
  306.  
  307. /****************************************************************
  308.  
  309.     or_init()        Set up an orbit
  310.  
  311. ****************************************************************/
  312.  
  313. or_init( orbit, ap, aa )
  314.     int orbit;
  315.     double aa, ap;
  316.     {
  317.  
  318.     /***    Calculate semimajor axis        */
  319.  
  320.     sfs_orbits[ orbit ]->semimajor = ( aa + ap + ( 2 * fo_radius )) / 2;
  321.  
  322.     /***    Calculate orbital center - focus distance       */
  323.  
  324.     sfs_orbits[ orbit ]->dist = sfs_orbits[ orbit ]->semimajor - ( fo_radius + ap );
  325.  
  326.     /***    Calculate semiminor axis        */
  327.  
  328.     sfs_orbits[ orbit ]->semiminor = sqrt( (double) ( sfs_orbits[ orbit ]->semimajor * sfs_orbits[ orbit ]->semimajor ) - (double) ( sfs_orbits[ orbit ]->dist * sfs_orbits[ orbit ]->dist ));
  329.  
  330.     /***    Calculate orbital eccentricity  */
  331.  
  332.     sfs_orbits[ orbit ]->eccentricity = sqrt( 1.0 - (( sfs_orbits[ orbit ]->semiminor / (double) sfs_orbits[ orbit ]->semimajor ) * ( sfs_orbits[ orbit ]->semiminor / (double) sfs_orbits[ orbit ]->semimajor )) );
  333.  
  334.     /***    Calculate orbital period        */
  335.  
  336.     sfs_orbits[ orbit ]->period = sqrt ( ( ( 4.0 * ( PI * (double) PI ) ) / ( UGC * fo_mass ) )
  337.            * ( sfs_orbits[ orbit ]->semimajor * sfs_orbits[ orbit ]->semimajor * sfs_orbits[ orbit ]->semimajor ) );
  338.  
  339.     /***    Calculate the increment factor of longitude of 
  340.         the ascending node.  This factor must be multiplied
  341.         by a time factor (orbital period or time into orbit, 
  342.         in seconds.  See Davidoff, pp. 8-9 and 8-10, and 
  343.         formula 8.14.     */
  344.  
  345.     if ( fo_siday == 0.0 )
  346.         {
  347.         sfs_orbits[ orbit ]->lif = 0;
  348.         }
  349.     else
  350.         {
  351.         sfs_orbits[ orbit ]->lif = (2 * PI) / (double) fo_siday;
  352.         }
  353.  
  354.     }
  355.  
  356. /****************************************************************
  357.  
  358.     or_kep()        Solve Kepler's equation
  359.  
  360.     Globals utilized:  
  361.  
  362.     sfs_orbits[ orbit ]->eccentricity Eccentricity of the orbital ellipse
  363.  
  364.     sfs_orbits[ orbit ]->semimajor    Semimajor axis of the orbital ellipse (km)
  365.  
  366.     sfs_orbits[ orbit ]->period       Orbital period (seconds)
  367.  
  368.     Inputs:
  369.  
  370.     t               Time from periapsis in seconds ( 0 < t < sfs_orbits[ orbit ]->period )
  371.  
  372.     Output:
  373.  
  374.     theta           Angle between periapsis, geocenter, and
  375.             current position (theta).  
  376.  
  377.     r               Distance from satellite to focal center,
  378.             in kilometers
  379.  
  380. ****************************************************************/
  381.  
  382. or_kep( orbit, t, theta, r )
  383.     int orbit;
  384.     long t;
  385.     double *theta;
  386.     long *r;
  387.     {
  388.     double z, z3, E;
  389.     z  = 2.0 * PI * ( (double) t / (double) sfs_orbits[ orbit ]->period );
  390.     E  = z;
  391.  
  392.     do
  393.         {
  394.         z3 = ( E - ( sfs_orbits[ orbit ]->eccentricity * sin( E )) - z ) / ( 1 - ( sfs_orbits[ orbit ]->eccentricity * cos( E )));
  395.         E  = E - z3;
  396.         }
  397.     while ( fabs( z3 ) > 0.000000001 );
  398.  
  399.     *theta = PI;
  400.     if ( E != PI )
  401.         {
  402.         *theta = 2.0 * atan( sqrt( ( 1 - sfs_orbits[ orbit ]->eccentricity ) / ( 1 + sfs_orbits[ orbit ]->eccentricity )) * tan( E / 2.0 ));
  403.         if ( E > PI )
  404.             {
  405.             *theta = ( (double) 2.0 * PI ) + *theta;
  406.             }
  407.         }
  408.  
  409.     *r = sfs_orbits[ orbit ]->semimajor * ( 1 - sfs_orbits[ orbit ]->eccentricity * sfs_orbits[ orbit ]->eccentricity ) / ( 1 + sfs_orbits[ orbit ]->eccentricity * cos( *theta ));
  410.     return 1;
  411.     }
  412.  
  413.  
  414. /****************************************************************
  415.  
  416.     or_ssp()    Calculate Subsatellite Point
  417.  
  418.     This function utilizes available data on the satellite 
  419.     to return the subsatellite point, that is, the point 
  420.     on the surface of the orbital focus directly under 
  421.     the satellite at a given moment.  
  422.  
  423. ****************************************************************/
  424.  
  425. or_ssp( orbit, t, latitude, longitude, r, n )
  426.     int orbit;
  427.     long t, *r, *n;
  428.     double *longitude, *latitude;
  429.     {
  430.     static long _r, tp;
  431.     static double theta;
  432.     double n1, n2, n4, E, lan;
  433.     long t1;
  434.  
  435.     *n = t / sfs_orbits[ orbit ]->period;                  /* return number of current orbit */
  436.     t1 = t - ( sfs_orbits[ orbit ]->period * ( *n ) );     /* get seconds into this orbit */
  437.     if ( t1 == 0 )
  438.         {
  439.         t1 = 1;
  440.         }
  441.  
  442.     or_kep( orbit, t1, &theta, &_r );
  443.     *r = _r;
  444.  
  445.     /***    Calculate the longitude of the ascending node at 
  446.         the beginning of this orbit.  See Davidoff, p. 8-9
  447.         and 8-10, and equations 8.13a, 8.13b, and 8.14.    */
  448.  
  449.     lan = sfs_orbits[ orbit ]->lon_an - ( sfs_orbits[ orbit ]->lif * sfs_orbits[ orbit ]->period * *n );
  450.  
  451.     /***    Calculate the latitude of the SSP.  See Davidoff, 
  452.         pp. 8-13 - 8-15, esp. equation 8.20.        */
  453.  
  454.     *latitude = asin( sin( sfs_orbits[ orbit ]->inclination ) * sin( theta + sfs_orbits[ orbit ]->arg_per ));
  455.  
  456.     if ( ( sfs_orbits[ orbit ]->inclination * RAD_DEG >= 0 ) && ( sfs_orbits[ orbit ]->inclination * RAD_DEG < 90 ) )
  457.         {
  458.         n2 = 1;
  459.         }
  460.     else
  461.         {
  462.         n2 = 0;
  463.         }
  464.     if ( ( ( *latitude ) * RAD_DEG ) < 0.0 )
  465.         {
  466.         n4 = 1;
  467.         }
  468.     else
  469.         {
  470.         n4 = 0;
  471.         }
  472.  
  473.     if ( ( sfs_orbits[ orbit ]->arg_per * RAD_DEG > 180 ) && ( sfs_orbits[ orbit ]->arg_per * RAD_DEG <= 540 ) )
  474.         {
  475.         n1 = 1;
  476.         }
  477.     else
  478.         {
  479.         n1 = 0;
  480.         }
  481.  
  482.     E = 2 * atan( (double) pow( ( 1 - sfs_orbits[ orbit ]->eccentricity ) / ( 1 + sfs_orbits[ orbit ]->eccentricity ), (double) 0.5 )
  483.               * tan( sfs_orbits[ orbit ]->arg_per / 2.0 )) + ( 2.0 * PI * n1);
  484.     tp = ( sfs_orbits[ orbit ]->period / ( 2.0 * PI ) ) * ( E - sin( E ) );
  485.  
  486.     *longitude = or_lfix( lan 
  487.              - pow( (double) -1.0, n2 + n4 )
  488.              * acos( cos( theta + sfs_orbits[ orbit ]->arg_per ) / cos( *latitude ) )
  489.              - sfs_orbits[ orbit ]->lif * (double) t1
  490.              - sfs_orbits[ orbit ]->lif * (double) tp );
  491.     }
  492.  
  493. double
  494. or_lfix( l )
  495.     double l;
  496.     {
  497.     double l1;
  498.     l1 = l;
  499.     while ( l1 < ( 0 - PI ) )
  500.         {
  501.         l1 += ( 2.0 * PI );
  502.         }
  503.     while ( l1 > PI )
  504.         {
  505.         l1 -= ( 2.0 * PI );
  506.         }
  507.     return l1;
  508.     }
  509.  
  510. ------------------------------
  511.  
  512. End of SPACE Digest V10 #89
  513. *******************
  514.